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The partition function of QCD is analyzed for an arbitrary number of flavors, Nf , 
and arbitrary quark masses including the contributions from all topological sectors 
in the Leutwyler-Smilga regime. For given Nf and arbitrary vacuum angle, 9, the 
partition function can be reduced to Nf — 2 angular integrations of single Bessel 
functions. For two and three flavors, the 9 dependence of the QCD vacuum is studied 



' in detail. For Nf — 2 and 3, the chiral condensate decreases monotonically as 9 

increases from zero to tt and the chiral condensate develops a cusp at 9 = tt for 
degenerate quark masses in the macroscopic limit. We find a discontinuity at 9 = tt in 
the first derivative of the energy density with respect to 9 for degenerate quark masses. 
This corresponds to the first-order phase transition in which CP is spontaneously 
broken, known as Dashen's phenomena. 

^ ■ I. INTRODUCTION 

> : 

For Nf flavors of massless quarks, the Lagrangian of quantum chrornodynamics (QCD) is invariant 
under the global symmetry group SU (Nf) ^ x SU(Nf) l X ^7(1)^4 at the classical level. At the quantum 
level, however, the SU (Nf)n x SU(Nf) l symmetry is spontaneously broken to the diagonal subgroup 
of vector transformations, SU(Nf)R+L = SU(Nf)v , by an expectation value for the quark-antiquark 
condensate. The U(1)a symmetry is explicitly broken to Z(Nf)A by a non-vanishing topological 
susceptibility ffl. For scales well below A, the typical hadronic mass scale, the effects of the explicit 
breaking of the U(1)a symmetry can be essentially ignored and the dynamics of the theory are 
dominated by the Nf — 1 Goldstone bosons that arise from the spontaneous breaking of the chiral 
symmetry. For small, nonzero quark masses, these excitations become pseudo-Goldstone bosons. 

In addition to the quark masses, the QCD Lagrangian depends on another parameter, 9, the 
so-called vacuum angle. Though this parameter explicitly violates the discrete CP symmetry for all 
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non-integer multiples of tt, experimental results for the value of the neutron dipole moment constrain 
9 to be zero within a deviation of less than |^||- The reason why 9 = 0, however, is still poorly 
understood, and there is considerable theoretical interest in the physics of QCD at 9 ^ 0. The most 
well-known example is the Veneziano-Witten formula which in the large- N c limit relates the mass 
of the if meson to the second derivative of the pure Yang-Mills vacuum energy with respect to 9 
Another phenomenologically motivated example is the work of Refs. |6|-|j| in which it is argued 
that metastable states in which 9 is effectively nonzero may be created in ultrarelativistic heavy ion 
collisions. 

Physics at 9 ^= is inherently nonperturbative and one must rely either on effective theories 
or perform lattice gauge theory simulations E(3] . The latter approach is prohibitively difficult at 
present because at nonzero values of 9 the action is complex and standard importance sampling 
methods that are usually employed are no longer applicable. This difficulty is similar to that for 
nonzero values of the baryonic chemical potential, but there has been recent progress in solving 
the complex action problem in simpler models ]Il]Jl2]] . While many observables associated with the 
quark masses have been studied extensively in lattice gauge theory computations, properties of the 
QCD vacuum associated with 9 such as the topological susceptibility have only recently been given 
extensive attention by the lattice community ]T^-|20[|. So, most of our knowledge of the physics 
of QCD at 9 has been gleaned from effective theories. One of the most striking examples 
of novel physics is the emergence of two CP violating degenerate vacua separated by a potential 
energy barrier at 9 = tt known as Dashen's phenomenon pi] ]. While QCD is seemingly invariant 
under the discrete CP group for 9 equal to integer values of tt, this invariance is spontaneously 
broken in a first-order phase transition at 9 = n. Chiral effective Lagrangians in the large- N c limit 
were used to investigate this and other properties of the 9- vacua in Refs. J22j and [p3[ . Subsequent 
elaborations and refinements using this approach were made in Refs. [^4-29|. Additionally, random 
matrix models [ jBOpl] ] as well as numerical simulations of CP N ~ X models p3-pj5| have also been 
employed to investigate the physics of nonzero values of 9. 

In this work, we investigate the properties of QCD at 9 ^ in the Leutwyler-Smilga finite volume 
regime p6| , also known as the mesoscopic regime. Here, chiral perturbation theory is valid but the 
volume is taken such that the Goldstone modes are constant and their kinetic term can be ignored. 
This approach differs from those of Refs. [^2 - 29 1 , in which the large- N c is taken and the Lagrangian is 
studied to lowest order in the chiral fields. The Leutwyler-Smilga regime, on the other hand, enables 
exact, analytical predictions to be made. Indeed, this limit is realized in many lattice simulations 
p7j. Many of the low-energy aspects of QCD have been studied using effective chiral Lagrangians 
in a finite volume [[38]]. In Ref. Leutwyler and Smilga demonstrated that such theories contain 
information beyond pion dynamics and the spontaneous breaking of chiral symmetry, including the 
low-lying spectrum of the QCD Dirac operator and the topology of the gauge fields. The purpose 
of this work is twofold. First, we reduce the QCD partition function for Nf > 3 in the Leutwyler- 
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Smilga scaling regime including the contributions from all topological sectors to Nf — 2 angular 
integrations over single Bessel functions. The latter half of this work is devoted to examining the 
vacuum properties of QCD in the Leutwyler-Smilga scaling regime at nonzero values of 9. 

The partition functions for one and two flavors of quarks including contributions from all topolog- 



ical sectors were computed in Refs. ]39(] and ]40[. As emphasized in Ref. 41 1, the contributions to 
the partition function from sectors with nonzero topological charge are crucial in the finite volume 
regime. For example, if one considers the partition function at any fixed topological charge (with 
the exception of v = 0), then chiral condensate diverges as one of the quark masses approaches 
zero. This completely unphysical behavior is remedied after summing over all topological charges, a 
procedure which becomes significantly more complicated for three or more quark flavors. We show 
in the following that the summation over all topological sectors can be performed for arbitrary Nf 
and 9 with the reduction of the partition function to Nf — 2 angular integrations over single Bessel 
functions. 

This paper is organized as follows. In Sec. ||, we briefly review the reduction of the QCD partition 



function to a finite dimensional group integral in the Leutwyler-Smilga scaling regime. In Sec. HI 



we review the summation over all topological charges for the partition function for Nf =1 and 2. We 
then calculate the summation for arbitrary Nf and quark masses. While the full partition function 
for Nf — 2 has been known for some time |]3rj] , there has not been a detailed study of the vacuum 



properties of this theory for 9^0. In Sec. IV, we examine the 9 dependence of the chiral condensate, 
chiral susceptibility, topological density and topological susceptibility for two flavors. In Sec. [V], we 
extend this discussion to three flavors. Finally, we give concluding remarks in Sec. (v|. 



II. PARTITION FUNCTIONS IN THE LEUTWYLER-SMILGA REGIME 

Consider a SU (N c ) gauge theory with Nf flavors of fermions in the fundamental representation 
on a four-dimensional torus of volume V — L . The Lagrangian is 

C = ^2 -^ f (m-m f )^-i9 f; v f« v , (i) 
9 /=1 

where F® v = ^e M1/a ^F°^. The last term in eq. ([!]) is a total derivative and so does not affect the 
field equations or any of the perturbative aspects of the theory. The integral of this term over the 
four-volume is quantized and is given by 

V = 32^ / d ' X e Z • ( 2 ) 

This quantity, known as the topological charge, contributes a phase factor to the path integral and 
is associated with transitions between topologically nontrivial gauge field configurations. The full 
partition function is given by a weighted sum over the partition function for each topological sector: 
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Z^(9 : {m k })= ^ 6 zl Nl \{m k }) . (3) 

V — — oo 

The individual contributions are 

zi Nf \{m k }) = J [dA] u ]ldet[ip[A]- mf ]exp(- J^xCymIA]^ , (4) 

where the integration is taken only over gauge field configurations with topological charge v. 

In principle, all of the observables of QCD can be calculated by evaluating eq. (||) and its various 
functional derivatives. Typically, however, this can only be done in some approximation since the 
integration is taken over an infinite dimensional functional space and there are ultraviolet and 
infrared divergences. Formulating the theory on a discrete, Euclidean lattice has the advantage of 
eliminating the ultraviolet divergences on account of a nonzero lattice spacing, a, and the infrared 
divergences on account of a necessarily finite volume. A second approach is to analyze the partition 
function for the effective theory which describes the low-energy behavior. For QCD, this is chiral 
perturbation theory. Restricting the Euclidean four- volume, V — L 4 , to the range 

11 

- « L « , (5) 

A ITItt 

where A is the typical hadronic scale and m w is the mass of the Goldstone excitations, results in 
tremendous simplifications and exact analytical predictions are possible p8| , |36| ]. This is possible 
since the lower limit ensures that the partition function is dominated by Goldstone modes and the 
upper limit ensures that these modes are constant, i.e. that the kinetic term in the partition function 
factorizes from the mass dependent term. 

There have been a number of advances in recent years towards evaluating the QCD finite volume 
partition function. With the four-volume taken according to eq. (||), the partition function was 
shown in Refs. |38|]36f | to reduce to the finite dimensional group integration 

Z ( - N f\0,{m i },V)= [ dU exp (e l6/N fTrMuA] , (6) 

where M. = diag(mi, . . . , TO/v*) is the mass matrix for the quark fields and £ is the chiral condensate 
in the chiral and infinite volume limit. We shall henceforth refer to £ as the macroscopic chiral 
condensate. The integration is taken over the group manifold of the Goldstone modes, SU(Nf). 
Note that the dependence on the volume, the macroscopic chiral condensate and the quark masses 
is only through the dimensionless scaling variable 

[ii = T*Vmi. (7) 

The Fourier coefficients conjugate to Z^ N f\0, {/^}) are obtained by taking the Fourier transfor- 
mation of eq. (||), 



4 



2tt 



(8) 



They have been computed for arbitrary Nf and N c > 3 ■ 

2v (|M, V) - A({A1 2 }) 



and can be expressed as 



(9) 



where 



is a AT/ x iVf matrix and /„ is a modified Bessel function of order n. The denominator, 



A({/i?})=n^-/# 



(10) 



1>J 



is the Vandermonde determinant. For equal masses, /i, = Z v (\i) can be further simplified using 
properties of Bessel functions and determinants to give 



I v -i(jjl) I„(ji) 



(ii) 



Combinatorial formulas for generating the character expansions of the U(N) group were recently 
given in Ref. fl43| and used to very efficiently derive eq. (|9|) . The close connection to random matrix 
theory was firmly established in Refs. Q-Q. 

The vacuum properties of the theory are defined by eq. (||) and its derivatives. Recall that the 
first derivative of the logarithm of the partition function with respect to a parameter determines 
the mean value of the conjugate variable to the parameter, while the second derivative serves as a 
measure for the fluctuations around the average value. We focus on the derivatives with respect to 
the quark masses, m„ and the vacuum angle, 9. The chiral condensate for the i-th quark species is 

1 d 



z\^(6,{m k },V) 



with the average value 



VN f dm, 



N f 



\ogZ^(9,{m k },V) , 



EW>(0, {m fe }, V) = J2 (°> {™ k }, V) 



(12) 



(13) 



i=i 



This quantity should approach macroscopic chiral condensate, S, at 9 = and for large values of 
the scaling variable, \i. The chiral (scalar) susceptibility is defined as 
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X<f' \e,{m k },V) = l^-l 0gZ W(e,{m k },V) . (14) 

Vacuum properties associated with the vacuum angle, 9, can be analogously defined. The topo- 
logical density is 

a(9,{m k },V) = -^logZ^(8,{m k },V) , (15) 

which at 9 — has the interpretation of the mean topological charge. From eqs. (^) and (^|), one can 
see that the topological density vanishes if 9 is an integer multiple of ir. The topological susceptibility 
is defined by 

Xu v (e ) {m t hV) = ~-^]ogZ^f\9,{m i hV) , (16) 
which at 9 = is the mean square deviation of the topological charge and is in general nonzero. 



III. SUMMATION OVER TOPOLOGICAL CHARGES 



When studying the vacuum properties, e.g. the chiral condensate, in the Leutwyler-Smilga regime, 
the necessity of including the contributions to the partition function from every topological sector 
was established in Ref. |4l}| . For one and two quark flavors, the summation can be performed using 
standard Bessel function identities. For Nf > 3, however, the summation is more complicated and 
to the best of our knowledge has not been performed previously in the literature. It is instructive to 
review the derivation of the results for one and two flavors before proceeding to the derivation for 
N f >Z. 



A. N f = 1 



For Nf — 1, the summation over all topological sectors is straightforward to calculate. Using the 
definition of the generating function for modified Bessel functions p7|. 



exp 



U— — QO 



the summation can be performed: 

oo 



(17) 



(18) 



1/— — 00 



exp [| (e^ + e- 9 ) 
exp [jj, cos 9] . 



This result was first found in Ref. |36[ 



G 



B. N f = 2 



For two flavors, the calculation is more involved but the method generalizes very naturally to 
Nf > 3. From eq. (|9|), the partition function is 



Z (N J= 2), 



1 — — 9 2 

M2 - Mi 



(19) 



f = — 00 



Since the partition function is symmetric in the two scaling variables, we consider only the first term 
between the brackets in eq. (n9J). Following Ref. p7|, we use the contour integral representation for 



I v (x) which follows from eq. (17), 



Iu{x) = <f> -^s 
' 2m 



exp 



2^ S + S 



(20) 



where the contour is the standard Bessel function contour surrounding the negative real axis and the 
origin in a counterclockwise fashion. Using this representation, the summation can be calculated: 



OC OO Jt j 



v= — co 

ds 
2ttI 

ds 
2ttI 



s exp 



s exp 



/'2 



-^(Mi) 



(s + s x ) I exp I ^— I se 



[Mi 

It 



s" 1 e ie 



By making the change of variables, 



(M2 +Mie~ 49 ) 
w = s 

M12W 
_i (M2 +/iie lS ) _ x 

where 

^12(6*) = \/mi + M2 + 2 Mi M2 cos 
eq. (|21a[) can then be simplified to yield 



£ e"^(MiK + i(M2) = 



M2 + Mie" 



Mi2W — .f 2™ 

-i6 



uj exp 



M12W 



(21a) 
(21b) 
(21c) 

(22a) 
(22b) 

(22c) 
(23) 



M2 + Mi e 



■Jl 0*13(0)) • 



M12W 

Combining the two contributions, the partition function for Nf = 2 is 

cr(Nf=Z)ra ^ MMi2(0)) 

M12W 



(24) 



Equation (22c) defines a reduced mass determined by a triangle law. The triangle has sides fix, ^ 2 
and /X21 and the angle subtended by \i\ and (j,2 is tt — 0. 
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C. Nf > 3 



For Nf > 3, we consider cases of degenerate and nondegenerate quark masses separately. In the 
case of nondegenerate quark masses, the summation over topological charge is a generalization of the 
summation for Nf = 2. In the process of the derivation, we find interesting relationships between 
sums of products of arbitrarily many Bessel functions and angular integrations over single Bessel 
functions. 

1. Nondegenerate quark masses 



From eq. (j|) , the quantity that needs to be calculated is a summation of products of Nf modified 
Bessel functions weighted by a phase, since 



detX({M) =£ii...i Nf '[[t4~ 1 I v +j-i(lM j ) ■ 
Each term in the summation is of the form 

oo 

6= e we /„ +mi (a;i)/, /+m2 (a;2) • ■ ■I v+mNf (x Nf ) 



(25) 



(26) 



J2 e we I v { Xl )h +n2 {x 2 )...I u+nN} {x N} ) 



U— — OC 



-irriiO 



c 



where rij = m,j — m,\. For the calculation of the partition function, B is actually more general than 
necessary and the particular choices of Xj = fiy and rrij = j — 1 specialize to the appropriate terms 
in the partition function. 

By expressing all but the first of the modified Bessel functions by their contour integral represen- 
tations, C can be written as 



C= £ e^^n 



dsj -v- 



n-j — l 



(27) 



The definition of the generating function for I n (x), eq. (|17[), can be used to perform the summation 
over v. 



N f 
3=2* 



dSj —n*j — 1 

-^s, exp 
2m J 



exp 



s 2 ... s Nf 



2 \s 2 ---s Nf 



We now make a change of variables which is suggested by eqs. (|22a|) and ( |22b[ ) 

S3 . . . Sn } \ S2 



U = \X2+ Xi- 



LJ 1 = | X 2 + X\ 



e" 



S 3 . . . SNf J 1p 



(28) 

(29) 
(30) 
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where 



'4> 2 = X x + x\ + X\X 2 



s 3 . . . S Nf 



s 3 . . . S Nf 



The contour integration over s 2 is rewritten as 



dS2 -n 2 -l 



exp 



2m A 

I x 2 + xis 3 . . . SN f e 



s 2 ( s 3 . .. S Nf \ 1 / < 
— [x 2 + xi — e ) + — ( x 2 + Xi 



2s 2 



S3 . . . SNf 



(31) 



(32) 



■n \ "2 



x 2 + X1S3 . . . SN f e~ 



The expression for C is then 

N 

r ( ds J 

2m" 3 



3=3 



x 2 + xis 3 . . . s^e" 



4 2 «0- (33) 



Making the second change of variables, Sk = e lc ^ k , for k = 3. . . , JVy, and deforming the contours to 
the unit circle, the final expression for C is obtained: 

"2 



2* 



i=3 



c = / n ( ^ exp ^ cos (^) - in 3$i\ 



\Jx 2 + x\ + 2x 2 xi cos((/> 3 + . . . + 4> 



N, 



x I n 



i2 I V -"2 



+ 2X2^1 cos( 



Collecting the various contributions, the partition function for Nf quark flavors is 



1 (r\ k -i\ 



where 



(34) 



and 



' 2;r 



-7T- exp HO - Wfc + Mi* cos(0 fe )] — — 



fc=3 



^Hi 2 (8,4>3, ■ ■ -,<f>N f ) = \J i"| + Mi 2 + 2//^/^ COS (<fo + . . . + (j) Nf - 1 



(35) 



(36) 



Note that S(fj.i 1 , /ii 2 ; Hi 3 , . . . , /Ltj^ ) is symmetric under the interchange of fj,^ and /ii 2 . 

In the case of complete nondegeneracy, fa ^ [ij for all i ^ j, the partition function and its 
derivatives are analytic functions to all orders. Any potential nonanalyticities in eq. p3) and its 
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derivatives occur when (1,^^(9; (f>3, ■ ■ ■ , 4>Nf ) vanishes which is only possible if at least two quark 
masses become degenerated. This becomes apparent when eq. (pq) is rewritten in the form 




which has a cusp at <fo + . . . + 0^ — = 7r. This behavior is a necessary condition for the derivatives 
of eq. (p4) to have discontinuities. On account of the integration, this is not a sufficient condition. 



2. Degenerate Quark Masses 

The limit of equal quark masses, [ii — > n, can be calculated either by summing over all topological 
sectors for nondegenerate quark masses and then taking the limit fii —> /i, or starting from the equal 
mass partition function in a given topological sector, eq. ([fl]) and then summing over all topological 
charges. Using the first approach, one can derive equal mass limit from cqs. (^)-(|36|) and the 
resulting expression for the partition function involves derivatives up to order Nf — 1 of modified 
Bessel functions with respect to the masses. This approach is quite cumbersome and, moreover, the 
derivatives are by construction taken at the cusps given by eq. (|37j). The second approach is more 
tractable and is the one used here. After expanding the determinant in eq. (pd|), 

Nf 

3=1 

and using the contour integral representation for each Bessel function, the expression for the equal 
mass partition function in a sector of given topological charge becomes 

Nf 

zi Nf \») = I II ex ? [ffo + *t)\ s j ' ' ' ) ■ ■ ■ 4? (39) 

3=1 

= / II (^T exp [mcos(^) - ii/^ - J e^n . . . e lN f^, . 

o i =1 

The last equality is obtained by change of variables, Sj — e l ^ j , deforming the integration contours 
to the unit circle, and making use of the identity 

e kl ... kNf e 1 ^ ...e N '^i = s kl ... kN / k ^ . . . e ikN f*»f . (40) 



zl N '\ii) = e h ... iNf J] lu+j-i, (P) , (38) 



^This is true for finite /ii. As shown for Nf = 2 in Ref. |4q] , there can be a nonanalyticity if the two 

scaling variables are taken to infinity such that lim^^^oo ^ ^ > We expect that this limit can be 

generalized to arbitrary Nf in a straightforward manner. 
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The authors of Rcf. [pw showed using a result of Weyl that 



Z " Nf)(p)= Nji.J a ll^exp[ M cos(^)+^i]Jll|e^-e^f . (41) 
Equations ( p9| ) and (|4l]) can be transformed into one another by observing that 

iNf<fr 

k<l 

is a Vandermonde determinant, and so 

Y[ |e l0fc -e #1 | 2 = A({e^})A({e-^}) . (43) 

k<l 

By expanding A({e" 1 ^}), one can transform both expressions into one another by an appropriate 
relabeling of the integration variables. 

The summation over the topological charges can be performed using 

/ ^p/to,...,^,,*) e* {9 ~* 1 -"~** t) = ffa + --. + 4>N f -9,<h,...,<l> Nf ,6), (44) 

provided / is 27r-periodic in <pi. For the subsequent integration over <^> 2 , an integral representation 
of the Bessel functions is used. The details of the calculation are given in the appendix. The result 
is 

2W(^)= £o ,. w f (a 2 - ai ) ft (P-^e^A JBl - B 'l^ 3 '"- , ^' g)) (45) 
x cos {(3 - ai - a 2 )(8 + <p 3 + . . . + <f>N f )/2 - (3 - 03)^3 - ... - (N f - a Nf )cj) Nf } , 



,2* W / 



|2 



EiL-Js,*** 1 ■ ■ ■ e N,4 " Nf = A({e^}) = J] (e** - e*') (42) 



where 

^(03,-.-,0w / ,6 | ) =M\/2 + 2cos(^3 + .-- + 0jv / -0) ■ (46) 

The summation over the completely antisymmetric tensor can be simplified to give a (Nf — 2)-fold 
integration over a sum of Nf — 1 Bessel functions multiplied by some phases, 

N f _ , 

jj. cos(0 m )— i(m— 3/2) 0n 

l 2tt ' 

m— 3 

r(-l) p 7- , ' o» r '(fo,...AN f ,0) ■ 47 



The phases Or ''fei ■ • ■ , , 0) are given by 



JV/— r 

n ^,...,^,0) = £ ^°' +r/2)(03+ --- + ^^^t(03,...,M), (48) 



S N fh 
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with 



/Jg^ (&,..., tor,) 



I ... e i(k-2)tj> 3 



a ik(f> 3 



i(N f -l)4> 3 



■y ... e i ( fc - 2 )0iv / e ikif> Nf e i{l-2)4> Nf il<pN f ... e »(JV/" 



(49) 



Despite appearances, this last expression is real. As mentioned above, the possible emergence of 
a cusp in the integrand occurs in the ratios of the Bessel functions and /^(</>3, . . . , <f>N f , &). At this 
stage, it is not a priori obvious that there are nonanalyticities after the integrations are performed. 
Thus, one must consider them on a case by case basis. This is done for Nf = 2 and 3 in the following 
sections. 



IV. VACUUM PROPERTIES FOR N F = 2 



Since the partition function for Nf = 2, eq. (|2J), can be expressed in closed form even after 
summing over all topological charges, calculating the vacuum observables is straightforward. Most of 
the work in the literature has focused upon the behavior of the chiral condensate and the topological 
susceptibility as functions of the scaling variable for either 8 = or in a sector of fixed topological 
charge. In this section, we demonstrate that QCD in the Leutwyler-Smilga scaling regime exhibits 
a great deal of nontrivial behavior at nonzero values of 9 for Nf = 2. 

A. Partition Function 

The partition function for two quark flavors is 



h (^1+^2+2^1^2 cos 9 j 

ZW= 2 )(^, M2 ) = \ ' ■ (50) 

y/flf + ^2 + 2 MlM2 COS 

At 9 = (2n + 1)tt, where n is an integer, it reduces to 

(9 = (2n + IK M 2 ) = h " M , 2|) , (51) 

i.e. the theory at 9 — 7r is equivalent to the theory at 9 = but taking one quark masses to be 
negative. This is in accordance with the standard lore. Taking the two quark masses to be equal, 
however, leads to the surprising result that the partition function is independent of the scaling 
variable, 

Z^= 2 \9 = (2n+l)n,^) = \ . (52) 
This can be seen also by starting from the equal mass partition function: 
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(m ,1, * h (uV2 + 2 cos 9) , 
/i V 2 + 2 cos 

Since Z^f^ 2 ^ (9, (j,i, /j, 2 ) is a smooth function of fj,±, \x 2 and 0, it does not matter in which order the 
limits of degenerate masses and 9 — (2n + l)ir are taken. This result may be understood in terms of 
the microscopic spectral density. As shown in Ref. Q], after summing over all topological charges, 
the full microscopic spectral density is related to a particular quotient of the partition function for 
Nf flavors at 9 and the partition function for Nf + 2 flavors at 9 + tt. Thus, the partition function for 
two flavors at 9 = it is related to the quenched (Nf — 0) theory at 9 = 0. This fact also demonstrates 
at a qualitative level that the partition function becomes independent of the quark masses if and 
only if Nf — 2 and 9 = it (49). It is apparent already at the level of the partition function that, at 
least in the Leutwyler-Smilga scaling regime, QCD is very different qualitatively at 9 — mr than at 
all other values of 9. 

This can also be understood directly from chiral perturbation theory. For Nf = 2, the leading 
order term in chiral perturbation theory, SRe [Tr {A4e l8 / Nf U^j] , vanishes at 9 = tt when M. = ml 
since the trace of an SU (2) matrix is always real. This term, however, is the only relevant term 
in the Leutwyler-Smilga scaling regime as can be seen from eq. (^). From this observation, it is 
clear that the partition function at 9 = tt is simply a constant in this case. The resolution to this 
problem was given in Ref. f24|j where higher order terms quadratic in masses were included in the 
effective chiral Lagrangian. A more extensive analysis in terms of chiral perturbation theory was 
given in Refs. |2^j29|| . By including higher order terms in the chiral expansion, the energy density 
then depends once again explicitly on the quark masses. These terms, however, are suppressed in 
the Leutwyler-Smilga regime and so we do not include them. 

The most fundamental quantity that can be derived from the partition function is the energy 
density which is defined by 

£(9,^,(12) = ~ log Z {N '- 2 \9^ 1 ^ 2 ) . (54) 

An expansion for large /ii2($) gives 

£{9^ 1 ^ 2 ) = -^r± + 0(\ogii l2 (9)) . (55) 

Note that the volume dependence drops out in the leading order term. In the macroscopic limit, 



-£mi-£TO 2 cos0 + (9(TO2/mi) , m! » m 2 
hm £{9, mi, ma, V) = \ (56) 

— 2mE| cos(0/2)| , mi — iri2 = m 



There is a constant shift, £q = -Emj, in the energy density which we subtract in the following. The 
energy density in both limits is plotted in Fig. [l]. It is minimized at 9 = and maximized at 9 = ir. 
For degenerate quark masses, there is cusp at 9 = n. 
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FIG. 1. The energy density for Nf — 2 in the macroscopic limit. The lower curve corresponds to degen- 
erate quark masses, while the upper curve corresponds to an infinite mass splitting. The region between the 
two curves contains all finite mass splittings. The global minimum of the energy density for all cases occurs 
at e = 0. 



B. Chiral Condensate 



The chiral condensate is calculated from eq. ( |l2] ) and for the first quark flavor is 

(N f =2) /a , S Ql +M2 COSfl) I 2 (//12(A)) 
2j (1) (0,//i,/Z 2 ) = — i-t; — • (57) 

{ ' 2 //i 2 (») ii(//i2(0)) 

The condensate for the second quark flavor is given by the same expression with /ii and // 2 inter- 
changed. The total chiral condensate given by the sum is 

v(W,=2)//) \ s (//1 +// 2 )(1 + cos6>) I 2 (//i 2 (6>)) 

For = 0, the individual contributions are equal. For 8 = n, however, they are equal in magnitude 
but opposite in sign: 

„(JV /=2 ), v . , . E7 2 (//i 2 (7r)) (JV /= 2), . 

(7T,Ml^2) =Slgn(//i -// 2 )-— ; r -r T = -S) 2 / 7T>/*1.M2), 59 

1 j 2 Jl(//l 2 (7r)) ^ 

with their sum, the total chiral condensate, vanishing for any value of the quark masses. 
Taking of the limit of degenerate quark masses, the chiral condensate is 

//) = E ^±^Z 72 (M^ + gg) (60) 
(1 < 2) V ' 4 / x ( M V2 + 2 cos 6>) V ' 
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Three observations may be made based upon eq. ( j60| ) about the physics of QCD in the Leutwyler- 
Smilga scaling regime at nonzero values of 9: 

• For fixed /i, the chiral condensate decreases monotonically in the interval 9 £E [0,7r) and in- 
creases monotonically for 9 € (ir,2ir). 

• As^ = mSy->co, the chiral condensate develops a cusp at 9 = it. 

• For any value of /i, the chiral condensate vanishes identically at 9 = tt. 

This behavior is shown in Figs. and |3|. The last observation follows of course from the constancy 
of the partition function for Nj =2 at 9 = n and is not expected to be true outside the Leutwyler- 
Smilga regime. 

1 
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FIG. 2. The chiral condensate in the macroscopic limit as a function of 9 normalized to the infinite volume 
value, E. The lower curve corresponds to asymptotically large quark mass splitting and the upper curve 
is limit of degenerate quark masses. The area between the two curves contains all finite mass splittings 
between the quark masses. There is a cusp at 9 = n for equal quark masses. 

For small values of the scaling variable, the chiral condensate has the expansion 

SC^KM = cofm) ^ _ cos^ ^ 3 + cc^/2) ^ + 0{[i7) 
Yl 2 12 48 

As expected, the chiral condensate is linear to lowest order in the quark mass and vanishes term-by- 
term at 9 = (2n + The expression for the chiral condensate simplifies in the limit of degenerate 
quark masses and in the limit of large quark mass splitting: 
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^ N i= 2 \9,m 1 ,m 2 ,V) = < 



— v2 + 2cos# , 

2 /i(m£VV2 + 2cos£ 



mi » m 2 



mi = m 2 



(62) 



The volume dependence of the chiral condensate in these two limits is easily calculated. In the 
limit of small volume, the chiral condensate vanishes in both limits of degenerate quark masses and 
large mass splittings. In the macroscopic limit, V 3> 1/Em, the chiral condensate is 



lim £ (Ar ' =2) (6»,mi,m 2 ,K) = 



£cos 2 (6»/2) + C(m 2 /mi) , mi » m 2 
£|cos(0/2)| , m 1 =m 2 - 



(63) 



and is independent of the quark masses in both cases. Any nonzero quark mass splitting, mi ^ m 2 , 
interpolates between these two cases as shown in Fig. ||. The value of the chiral condensate is 
remarkably insensitive to the magnitude of the quark mass splitting. 



L 

J i I i L 

5 10 

FIG. 3. The chiral condensate for Nf = 2 as a function of (j,. The upper, middle and lower curves 
correspond to 8 = 0, 7r/2 and n, respectively. As 8 increases from zero to ir, the condensate decreases and 
vanishes identically at 8 = n. 

In Fig. |^, the chiral condensate as a function of the scaling variable is plotted for 9 = 0, n/2, and 
7r. As 9 increases from zero to 7r, the condensate decreases monotonically and vanishes identically 
at 9 — 7T. The monotonic decrease in the chiral condensate can be understood from the underlying 
QCD Lagrangian fl^| . If one performs a chiral rotation by an angle a, then 9 — > 9 — a and the chiral 
condensate, (ipip) rotates into (^75^). For Nf = 2, the (scalar) chiral condensate is completely 
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rotated into the pseudoscalar chiral condensate at 9 = n. This can be demonstrated analytically 
by introducing a complex mass matrix as a source for both the scalar and pseudoscalar chiral 
condensates. 



C. Chiral Susceptibility 



The chiral susceptibility has been previously studied in Ref. Q but the analysis was restricted 
to fixed values of v. As shown in Ref. ^fj, however, the contributions from all topological sectors 
can be important in the mesoscopic regime. Defining 



i HmW) 



m{0)h{m{0)) 

the diagonal terms in the chiral susceptibility are 

(Mi + cos flV ' 



The off-diagonal terms are 

X12 = X21 = 



M12W 2 

(fl 2 + Mi cos6>) 2 



(Mi + H2 cos 9) (/x 2 + Mi cos 9) 



S{^{6))+r{ m {6)) 
S^ 12 {9))+T(pi 12 {9)) 



S(^ 12 (9))+cos9T(p 12 (9)) 



M12W 2 

In the limit of degenerate quark masses, mi = m 2 = m, chiral susceptibilities become 



X11 = X22 = - [ficos 2 (9/2)S(fi(9)) + iSr(n(0))] 
m 

S 

X12 = X21 = — 
m 



[H cos 2 (0/2)<S(/z(0)) + n cos 0T(n{0)j\ 



At 9 = tt, the diagonal and off-diagonal elements are maximally different, 



Xn(0 = tt) = X22(# = tt) = 



4m 



X12(A = TT) = X2l(<? = TT) = -— = — 



An expansion in small /j, yields 



X11 = X22 = — 
m 



X12 = X21 = — 
m 



4 24 

On the other hand, in the limit /i> 1 for 9 =/= n 



^-lcos 2 (0/2)(2 + cos^ 3 +O( M 5 ) 
^ C0S9 cos 2 (9/2) ( 1 + 2 cos 9)[i 3 + O (fi 5 ) 



(64a) 
(64b) 

(65a) 
(65b) 

(66) 

(67a) 
(67b) 

(68a) 
(68b) 

(69a) 
(69b) 
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X22 = — 

m 



3 cos 9 



Xl2 = X21 



E 
rn 



sin 2 (6>/2) 

2|cos(0/2)| ^ 8/^ cos(0/2) 2 
" sin 2 ((9/2) 3 



o(i/ M 2 ) 

0(1/ M 2 ) 



(70a) 
(70b) 



_2|cos(0/2)| 8^cos(0/2) 2 

therefore, to leading order in 1/fi, Xu = X22 = — X12 — — X2i- The chiral susceptibility was shown 
to be proportional to 3/8m 2 V at = in Ref. EG]. The importance of summing over all topological 



charges is clear in this expression, since the leading order terms in eqs. (|70a| ) and ( |70b| ) vanish at 
= 0. 



D. Topological Density: First— Order Phase Transition at 6 = n 

The topological density, defined by eq. (|l5|), is 

fa \ 1 MiM2sin0/ 2 (/ii 2 (0)) , . 

a(0, Ml , M2 ) = - Ji(/|u(fl)) . (71) 

At = 0, er(0, fi±, /x 2 ) has the strict interpretation as the mean topological charge. Since in a large 
enough ensemble of gauge fields the average number of instantons should be equal to average number 
of anti-instantons, the topological density should be zero at = 0, and, indeed, the right hand side 
of eq. ( |7l| ) vanishes for 9 equal to any integer multiple of n. 

In the limit of degenerate quark masses, the topological density becomes 



1 Msm0 

V V2 + 2cos6> h (fiy/ 2 + 2 cos 0) 

The most interesting property of this relation is that in the limit of very large scaling variable, 
fi 3> 1, cr(0, /1) develops a discontinuity at 9 = n. This is the first-order phase transition proposed 
by Dashen pi. 



(j(9,mi,m 2 ,V) = < 



Em 2 sing r 2 ^„ T , — ^ +C(m 2 /mi) , mi » m 2 
7i(LI/toi) 

^ . ,„ _I 2 (2EVmcos(0/2)) 

V 7 ; ii(2EVrocos 0/2 ) ' 1 2 



(73) 



Then, in the infinite volume limit, 



. Em 2 sin0 + C(m 2 /mi) , mi > m 2 
hm cr(0,mi,m 2 ,y) = { (74) 

Em sin(0/2)sign(cos(0/2)) , mi = m 2 = m . 



The limit mi 3> m 2 coincides with the one-flavor case. For two flavors of degenerate quarks, there 
is a first-order phase transition at 9 — n, while for two flavors of non-degenerate quarks the phase 
transition disappears as long as the quark masses are kept finite. As mentioned above, it was shown 
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in Ref. |4q ] that there exists a particular scaling limit in which there may be a phase transition for 
nondegenerate quark masses. 
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FIG. 4. The topological density of the 8 vacua. The solid line corresponds to the limit where one quark 
mass is very heavy, mi 3> m,2, the long dashed line to mass splitting of m\jm-i — 2, and the dashed line 
to the case of degenerate masses, mi = m2 = m. A first-order phase transition occurs at 9 = tt in the 
case of degenerate quark masses. The transition is washed out by mass splitting between the quarks. The 
topological density is measured in units of m2£, where m,2 denotes the mass of the lighter quark. 



E. Topological Susceptibility 

The topological susceptibility reads 

X {Nf = 2) (0,mi,m 2 ,V) (75) 

1/ 2 OM0)) MiMl ^ia(0)coB0 , ( 4 7 2 GM0)) M/M0)A • 2 / m 

sin (^(7 ) 



For equal quark masses, one has 

v ft= 2 ), ™ m = Sm 2-cosfl MKg)) , S 2 m 2 ^ sin 2 (g) / J 2 (/i(g)) 2 _ \ 

X [ > > V) 2 |cos(0/2)| I^)) + 4 cos2(^/2) U(M#)) 2 /' k ' 

This reduces at = n to 

X ^= 2 )(^^,m,y) = -^F, (77) 

which diverges in the infinite volume limit on account of the first-order phase transition. For 9 =/= n, 
however, the limit is finite, 
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lim x ( JV /= 2 ) (0, m , V) = cos(#/2) I , ^ tt . (78) 

This is completely consistent with the expectation from the flavor singlet Ward-Takahashi identity 
at 6 — which predicts a linear rise in topological susceptibility with mass J5l| , 

lim X (Nf \0 = 0, to, V) = ^ + 0{m 2 ) . (79) 



Equation ( |78|) may be considered to be a generalization of eq. ([79D to nonzero values of 6. 
In the limit of large mass splitting, the topological susceptibility reduces to 

X < - Nf=2 H0,rni,m2,V) = £to 2 cosfl ?\J^ r — 4 + 0(m 2 /m{) , mi > to 2 , 

7i(ST/mi) 

which in the macroscopic limit is 



(80) 



lim x {Nf ~ 2) ^,m 1 ,m 2 , V) = £ to 2 cosfl + 0(to 2 /toi) , to-i > m 2 . (81) 




V. VACUUM PROPERTIES FOR N F = 3 

The summation over all topological charges in the partition function for Nf > 3 was performed 
in Sec. [II. In this section, we complement the analysis of Sec. |W| by examining the 9 dependence 
of the QCD vacuum for three flavors of quarks. The partition function for Nf = 3 involves a single 
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integration which we were not able to calculate analytically. Since the partition function has a group- 
theoretic origin, we expect that the integrals that result after summing over all topological charges 
provide some functional representation of the Goldstone manifold. However, since the integrand is 
a smooth function of 9 and fi except possibly at a single point in the interval of integration and 
the integral is taken over a compact interval, the numerical evaluation of the partition function is 
straightforward. Additionally, the derivatives of the partition function with respect to 9 and /z can 
be commuted past the integration. Explicit expressions are only given for the chiral condensate and 
the topological density since the expressions for the chiral and topological susceptibilities become 
quite complicated and are not particularly enlightening even in the limit of degenerate quark masses. 
We focus on a triplet of quark masses, (mug^t, mughtj ^heavy), between the limits of total degeneracy, 
TOnght = m hc avy, and very large mass splitting, mu g ht < m hc avy 



A. Partition Function 



For three degenerate masses, the partition function can be written in the compact form: 

2tt 



z( N f= 3 \9,fi) = 



d<t> 



^ cos c, 





2tt 



d(i) 



MP, 



g/i cos <p 

m(M) 



cos 



/ 30- 6> 



2 



h{jJL{9,4>)) - cos 



30 



cos(^) cos(^) 2 



(82) 



where 



2/x 



■ Unlike the Nf = 2 partition function, we find that Z iN f= 3 ^(9,fi) 



does not become independent of /j, at 9 = n which is consistent with the discussion in Sec. IV A. As 
discussed in Refs. [ ^8|j29|l , the leading order term in the chiral expansion does not vanish at 9 = ir 
for N f >3. 

The dependence of the energy density as a function of 9, however, is similar to the two flavor 
case. The energy density in the macroscopic limit is shown in Fig. [| shifted by its zero point, 
£{) = — Emheavy The global minimum and maximum of the energy density in the macroscopic limit 
are at 9 = and tt, respectively. As long as two quark masses are degenerate, there seems to be 
cusp at 9 — 7T. In general, the positions for the global minimum and maximum are 9 = and 
7T, respectively, for any value of fx. The zero point energy is dominated by the heavy quark mass. 
After subtracting this physically irrelevant zero point contribution, the energy density is relatively 
insensitive to taking the heavy quark mass to infinity. 



B. Chiral Condensate 

For three degenerate quark masses, the chiral condensate is 
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7T Z( N I= 3 )(0,n) 



a fl cos <p 



axh{li{e,(t>)) + a2h{Li{6,<t>)) 



where 



«i 



a-2 



/3<j)-9 



cos m cos 



3</> -9 
2 

- 2 cos 



30- 



(83a) 

(83b) 
(83c) 



In many respects, the behavior of the equal mass chiral condensate for two and three quark flavors 
is again similar. The chiral condensate as a function of the scaling variable is plotted in Fig. 
The chiral condensate decreases monotonically as 9 increases from zero to tt like the Nf = 2 case. 
However, it is nonzero at 9 = tt unlike the Nf = 2 chiral condensate. For small values of /z, the 
chiral condensate is linear in fi and independent of 9. 
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FIG. 6. The energy density in the macroscopic limit as a function of 9. The upper (lower) curve corre- 
sponds to infinite (zero) mass splitting between the two degenerate quark masses and the third quark mass. 
All finite mass splittings are contained in the region between the two curves. 

Figure |^ shows the chiral condensate as a function of 9 in the macroscopic limit, /i — > oo. In this 
limit, Yl( N f~ 3 ' (6, /i) has a cusp at 9 = tt, but otherwise is a smooth function. In the macroscopic 
limit with two degenerate quark masses, the value of the chiral condensate is very insensitive to the 
mass of the third nondegenerate quark. For any nonzero value of 9, the chiral condensate is always 
less than the 9 = chiral condensate. We find numerically that 
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£(0 = 7T, miight, m h cavy) 



, TTT-light — ^hcavy 



- i ""Might "C TOhoavy 



(84) 



to a very high precision. 
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FIG. 7. The mass dependence of the chiral condensate for various values of 8 for three degenerate quark 
flavors. The upper, middle and lower curves correspond to 6 = 0, n/2, and ir, respectively. As 8 is increased 
from zero to it, the chiral condensate decreases in magnitude for a given \x. 

C. Topological Density: First-Order Phase Transition at 6 — n 



For degenerate quark masses, the topological density is given by 



tr(6,n,V) 



1 



1 



7T Z( N 1=V V 



(hp 



gfl COS <f> 

MM 



hh(fi(e,^) + b 2 i 2 ( f i(d,ct ) )) 



where the coefficients are given by 



b\ = cos 

&2 = COS 



V 2 



— sin 

/K COS 



2 

/3<£- i 



sm 



3 

— tan . 
2 V 2 



- tan [ ^— - 

2 V 2 



(85a) 

(85b) 
(85c) 



V 2 y ~~ V 2 

The topological density as a function of is shown in Fig. in the macroscopic limit. For three 
degenerate quark masses, there is a discontinuity at 9 — -k which is simply Dashen's phenomena. 
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This discontinuity is washed out for any nonzero mass splitting, however, for any mass splitting 
between mheavy and TOught, the transition is always extremely rapid. 
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FIG. 8. The three-flavor chiral condensate in the macroscopic limit as a function of 9 for which two of the 
quark masses are degenerate. The upper (lower) curve corresponds to infinite (zero) mass splitting between 
the two degenerate quark masses and the third quark mass. All finite mass splittings are in the region 
between the two curves. 



The topological susceptibility in the macroscopic limit for Nf — 3 is shown in Fig. |l0|. For three 
degenerate masses, the topological susceptibility diverges at 9 = n on account of the first-order phase 
transition. For three quark flavors with only two degenerate masses, there is still an extreme drop in 
the topological susceptibility at 9 = 7r even for a very large mass splitting. For two degenerate quark 
flavors, the topological susceptibility is positive except in the vicinity of 9 = n. We find numerically 
that 



D. Topological Susceptibility 



S mii g ht 



flight — ^hcavy 



3 



lim x { f 3 '(9 = 0, m^M, nihcavy, v ) = 



(86) 



2 



m-light <§C TOheavy 



which is consistent with the Ward-Takahashi identity (|79|). 
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FIG. 9. The Nf = 3 topological density as a function of 9 in the macroscopic limit for which at least 
two of the quark masses are degenerate. The solid, dashed and long-dashed curves correspond to the mass 
splittings mheavy S> wiiight, m-hcavy = 2miight, an d ?Tiheavy = miight, respectively. For completely degenerate 
masses, there is a first-order phase transition at 9 = n. Even for large mass splittings, however, there is still 
a very rapid crossover. 
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FIG. 10. The Nf — 3 topological susceptibility as a function of 9 in the macroscopic limit for which at least 
two of the quark masses are degenerate. The solid, dashed and long-dashed curves correspond to the mass 
splittings m heavy > mi ight , m hoavy = 2mi ight , and m hoavy = ranght, respectively. For completely degenerate 
masses, the topological susceptibility diverges at 6 = tt on account the first-order phase transition. 

VI. CONCLUSION 

In this paper, we investigated the properties of the QCD partition function in the Leutwyler- 
Smilga finite volume scaling regime. The full partition function including the contributions from 
all topological sectors for Nf — 2 has been previously calculated. We extended these results by 
performing the summation over all topological charges for arbitrary Nf. For Nf > 3, the partition 
function can be expressed as a (Nf — 2)-fold angular integration over a finite sum of modified 
Bessel functions. We considered both the cases of degenerate and nondegenerate quark masses. The 
partition function remains 27r-periodic in 6 after summing over all topological charges. 

We systematically investigated the 6 dependence of the QCD vacuum in the Leutwyler-Smilga 
regime. In this limit, the partition function only depends on 9 and the scaling variables, /Ltj = mj£V . 
In the limit of degenerate quark masses, the Nf = 2 partition function is independent of the scaling 
variable at 9 = tt. As a consequence, the chiral condensate vanishes identically at 9 = tt for all 
values of fi. For fixed jj,, the chiral condensate decreases monotonically as 9 is increased from to tt. 

In the macroscopic limit, i.e. for \i — > oo, the behavior of the two-flavor partition function is 
particularly striking. For degenerate quark masses, the first derivative of the energy density with 
respect to 9 has a discontinuity at 9 = tt corresponding to the spontaneous breaking of the discrete 
CP symmetry. This phenomena was predicted by Dashen and was subsequently demonstrated by Di 
Vecchia and Veneziano and Witten using large- N c chiral perturbation theory. The chiral condensate 
also develops a cusp at 9 = tt for degenerate quark masses in the macroscopic limit. 

For Nf = 3, we find that all the examined quantities are very insensitive to the mass splitting 
between two degenerate flavors and the third heavy flavor in the macroscopic limit. When compared 
to those of Nf = 2, the vacuum properties for three flavors are quite similar. While the partition 
function for Nf = 2 can be expressed in closed form, the partition function for Nf — 3 requires 
the straightforward numerical evaluation of a single integral. Unlike the Nf = 2 case, the Nf = 3 
partition function is not independent of the scaling variables at 9 = tt and subsequently the chiral 
condensate is nonzero at 9 = tt even in the macroscopic limit. The chiral condensate does, however, 
exhibit the same monotonicity in 9 as the Nf = 2 chiral condensate and in the macroscopic limit 
there is a cusp at 6 = n. Dashen's phenomena, a first-order phase transition at 9 = tt, is also realized 
in this limit. 

Many of our results are corroborated by general 9 = field theoretic identities, for example, 
the Ward-Takahashi identity. Our analysis naturally extends these results to nonzero values of 9. 
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We have also demonstrated that QCD in the mesoscopic regime exhibits Dashen's phenomena. We 
examined the physics of QCD at 6 ^ in a nonpcrturbative framework. While the existence of 
Dashen's phenomena has been known for thirty years and studied using chiral perturbation theory 
and numerical simulations, we have demonstrated the spontaneous breaking of the CP symmetry 
at 6 = 7r in a way that is both nonperturbative and analytic. We have also provided a further step 
towards the full evaluation of the path integral for chiral perturbation theory. We hope that this 
work helps elucidate a parameter space of QCD which has been largely unexplored. 
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APPENDIX A: EXPLICIT CALCULATION OF THE DEGENERATE MASS PARTITION 

FUNCTION 

Starting with the representation 

27r Nt 



zi N '\») = J n exp [mcos^) - ivfr .«.v/' 0i ' ■ ■ ■ , (Al) 

o i =1 

the summation 2W(«,/<) = ^e^Z^fji) can be performed by the use of 

/ SgLffa,...,^^) e* {e -* 1 "--+ r,/) = f(<h + --- + <t>N,-9,<h,...,4>N f ,e), (A2) 

v=— oo 

After a change of integration variables 

e i<t>2( 1 + e i(4>3+-+4>N f -e)^ = ^ e i<p ^ 

e -*( 1 + e -'(fc+-.+fe r «)) = ^ e -i V ( A4 ) 

with 



V> = ^2 + 2008 (03 + . . . + 4> Nf - 9) (A5) 



one arrives at 



2^(6,1*) (A6) 

2-k 2ir 

= e -i9/2 J d^^cop-up J JJ^eMcoBWiMCi-a/^^r^^...^^^) 
n n J =3 
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where 



T(cp, <j> 3 , . . .,<j) Nf ,6) = 



I e ~H'l > 3---+<t>Nj,-0)-iifi 



i(N r —1) 

e - 2 - {4>3...+4>H f -e)-i{N J -i) v 



<j> 3 ...+<f> Nf -e)+i{N f -i) v 

e i(N f -l)<p 3 



J{N f -l)^ N 



Expanding the determinant around the first two rows leads to 



Nf-l 



a 



(Nf), 



i, ■ ■ ■ , <rN f i "i j 



where the phases ctr are given by 



Nf-r 



J'=l 



with 



4^ (03,...,^) 



1 ... g*(fe-2)03 g ife03 . . . gj('~ 2 )03 . . . gt(JV/ — 1) 



... e i(k-2)4> Nf e %k4>m t . . . g2('-2)0jv / e il<t>N f . . . e i(Nf-~L)<t> 

The integration over the angle </? can be easily performed, since 



2^' 



2r 

= rlr(^ip) ■ 

flip 



From that the final expression for the partition function with degenerate masses is 

2 7T Nf 



27T N f / J J. \ ^z -1 



<Wm ucos(<j. Ttl )-i(m-3/2)0 i 
m— 3 r— 1 



Alternatively, from cq. (11), the equal mass partition function is 

oo Nf 



Z^(9^)=e ai ... aNf E ^ II W-a m (M) 

1/— — OO 771=1 

= e 01 ...a JV/ E ^ II / ~~^T exp [MCOs(<?!) m ) + i(v + m - a 

v— — oc m— 1 . 



m/rm 
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Using eq. (p4|), the summation over v gives a delta function which facilitates the integration ai : 
Z( N f\0, n) = ^'jvT-i / d< ^ 2 • ■ ■ d< ^/ ex P [mcos(6» + 2 + 03 + ■ • • + <f>N f j\ (A14) 



x exp [/zcos(0 2 ) + M cos (03) + • • ./xcos(0jv>)] 

x cos{(ai - 1)(0 + 02 + . . . + 0^) + (2 - a 2 )0 2 + . . • + (N f - a Nf )cj) Nf } . 

The integration over 2 can be done analytically [p2| by expanding cos(# + 02 + . . . + 4>N f ) = 
cos(0 2 ) cos(6l + 03 + . . . + (j) Nf ) - sin(0 2 ) sin(0 + 3 + . . . + (f> Nf ). 
The result is 

n f J^e^A I aM (2 M cos( g + ^ + 2 - + ^ )) (A15) 
x cos {(3 - ai - a 2 )(6» + 03 + ... + (j>N f )/2 - (3 - a 3 )0 3 - ... - (JV/ - a Nf )<p Nf } . 
Using the addition rules for modified Bessel functions, one arrives at eq. (fl5|). 
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